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На основе разработанной математической модели проведен вычислительный эксперимент по исследованию динамики процес- 
сов теплопереноса в Н-образном информационно-логическом блоке современного космического аппарата связи при его орби- 
тальной эксплуатации в существенно нестационарных условиях суточного цикла освещенности Солнцем. Сделана количествен- 
ная оценка степени неоднородности температурных полей как приборов радиоэлектронной бортовой аппаратуры, так и основ- 
ных элементов конструкции приборного отсека. 


Современные искусственные спутники Земли, 
находящиеся на различных орбитах, используют 
для связи и телерадиовещания, дистанционного 
зондирования Земли, исследования ресурсов моря 
и суши, метеорологии и других целей. Они пред- 
ставляют собой долгоресурсные космические ап- 
параты (КА), в процессе эксплуатации находящи- 
еся в условиях комплексного воздействия внутрен- 
них и внешних тепловых нагрузок. Приборный от- 
сек (ПО) такого КА представляет собой сборную 
негерметичную бл очно -модульную конструкцию 
из прямоугольных трехслойных алюминиевых па- 
нелей с сотовым заполнителем [1] и собран из мо- 
дуля полезной нагрузки (МПН) и модуля служеб- 
ных систем. Модуль служебных систем состоит из 
вложенных друг в друга Н-образного информа- 
ционно-логического блока (ИЛ Б) и П-образного 
энергодвигательного блока (ЭДБ) [2], рис. 1. 


Внешний теплообмен ПО характеризуется пе- 
риодическим полугодовым циклом освещенности 
Солнцем (панели «Север» и «Юг»), Применение 
Н -образной конструкции ИЛ Б позволяет дополни- 
тельно использовать панели «Восток» - (1) и «За- 
пад» - (3) в качестве радиаторов. По этой причине 
внешний теплообмен в ИЛБ имеет существенно не- 
стационарный характер в суточном цикле поперемен- 
но-периодического воздействия солнечного потока. 

Ранее было проведено численное моделирова- 
ние пространственных нестационарных темпера- 
турных полей перспективных МПН и ЭДБ [3, 4], ра- 
ботающих в условиях стационарных внешних воз- 
действий. Для ИЛБ такой анализ проведен только в 
рамках плоской модели теплопереноса [5], которая 
не учитывает ряд важных особенностей Н-образной 
конструкции блока. Так, например, перенос тепло- 
вой энергии в такой системе по всем трем координа- 



Рис. 1 . Общая схема блочно-модульной компоновки перспективного КА связи 
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тным направлениям и связь между каждой из трех 
панелей играют важную роль в формировании те- 
плового режима блока. По этим причинам предста- 
вление каждой панели бесконечно тонкой пласти- 
ной или пластиной с бесконечно большой тепло- 
проводностью может привести к достаточно боль- 
шим погрешностям при численном моделировании. 


Цель работы заключается в математическом мо- 
делировании теплопереноса в информационно-ло- 
гическом блоке негерметичного ПО КА при его ор- 
битальной эксплуатации в существенно нестацио- 
нарных условиях суточного цикла освещенности 
Солнцем с учетом объемного характера переноса 
тепловой энергии. 

Задача рассматривается на примере перспектив- 
ного Н-образного ИЛБ, в котором приборной явля- 
ется панель «Центральная» - (2), имеющая двухсто- 
роннюю компоновку расположения приборов те- 
пловыделяющей бортовой аппаратуры (БА). Блок 
выполнен из алюминиевых трехслойных сотовых па- 
нелей (1-3) с встроенной системой обеспечения те- 
плового режима (СОТР) на основе нерегулируемых 
низкотемпературных тепловых труб (НТТ), связы- 
вающих в тепловом отношении панели - 1 и 3. Па- 
нели -1,3 открыты для теплообмена с космическим 
пространством. В качестве приборов рассмотрена 
типичная конструкция блока элементов радиоэлек- 
тронной аппаратуры кассетного типа с кондуктив- 
ными теплостоками, по которым осуществляется от- 
вод тепла от плат с микромодулями к охлаждаемому 
основанию. Мощности источников тепловыделения 
считаются равномерно распределенными по обла- 
стям, имеющим форму параллелепипедов и соответ- 
ствующим отдельным функциональным элементам 
или группам элементов [6]. В этом случае задача сво- 
дится к решению системы трехмерных уравнений те- 
плопроводности для N анизотропных параллелепи- 
педов с заданным объемным распределением вну- 
тренних источников теплоты: 


с,Р, 


дТ 
д I 


д_ 

дх 



ду 



+ 


д_ 
д г 



Ч, (х,у,г). 


Здесь с - удельная теплоемкость; р - плотность; 
Т— температура; 1 - время, Х„ А я А г - эффективные 
коэффициенты теплопроводности анизотропных 
элементов нагретой зоны; х, у, і - прямоугольные 
координаты; дІх,у,т) - функция, учитывающая ра- 
спределение мощности внутренних источников те- 
пла; /=1, 2 ,... ,Ц — индекс элемента нагретой зоны. 

Система уравнений теплопроводности для об- 
шивок панелей, сотового заполнителя и корпусов 
тепловых труб имеет аналогичный вид за исключе- 
нием источников тепла. На всех границах в местах 
контактов элементов внутренней структуры прибо- 
ра между собой и с охлаждаемым основанием, ос- 
нования прибора с обшивкой, обшивок с сотовым 
заполнителем приняты условия непрерывности те- 
пловых потоков и температур [3,4]: 
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Здесь А я А у , А, - эффективные коэффициенты 
теплопроводности анизотропного сотового запол- 
нителя. Для изотропных алюминиевых обшивок 
А,=А=А Г 

Радиаторные панели - 1 и 3 находятся в перио- 
дическом суточном цикле освещенности Солнцем. 
Можно считать, что космическое пространство по- 
глощает все уходящее с радиационной поверхности 
тепло как идеально черное тело с нулевой темпера- 
турой. Для внешних граней панелей - 1 и 3 прини- 
мались граничные условия, описывающие падаю- 
щий поток солнечного излучения в виде периоди- 
ческих аналитических зависимостей для панелей- 
радиаторов - 1 и 3 соответственно: 
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Здесь А 5 - коэффициент поглощения прямого 
солнечного излучения; 5 - плотность теплового 
потока солнечного излучения; <р - угол между нор- 
малью к поверхности панели-радиатора «Север» и 
«Юг» и направлением на Солнце; е - интегральная 
степень черноты поверхности; ст - постоянная 
Стефана-Больцмана; \ѵ - индекс поверхности па- 
нели-радиатора. 

Эффективный сброс избыточного тепла в окру- 
жающее космическое пространство осуществляет- 
ся за счет нанесения специального терморегули- 
рующего покрытия «оптический солнечный отра- 
жатель» на радиационные поверхности панелей 
(е»=0,85). Коэффициент поглощения прямого сол- 
нечного излучения изменяется во времени и зави- 
сит от степени деградации терморегулирующего 
покрытия при воздействии факторов космического 
пространства. Он изменяется от 0,08 в начале эк- 
сплуатации до 0,33 в конце 12-15-летнего срока 
активного существования КА [3]. Плотность те- 
плового потока прямого солнечного излучения 5 
зависит от положения Солнца при его видимом го- 
довом движении по большому кругу небесной сфе- 
ры (эклиптике). В точках зимнего и летнего 
солнцестояния 5=1440 и 1350 Вт/м 2 соответствен- 
но. В точках весеннего и осеннего равноденствия 
5=1395 Вт/м 2 [3]. 
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Математическое моделирование теплопереноса 
в нерегулируемых НТТ, встроенных в сотовый за- 
полнитель, проводилось с использованием запи- 
санной в цилиндрических координатах осесимме- 
тричной математической модели работы СОТР в 
распределенных параметрах, описывающей гидро- 
динамику и тепломассообмен в паровой и жидкой 
фазах рабочего вещества и кондуктивный теплопе- 
ренос в корпусе тепловой трубы [7]. 

На основе сформулированной математической 
модели был проведен вычислительный экспери- 
мент с целью исследования динамики процессов 
теплопереноса в Н-образном ИЛБ перспективного 
ПО современного КА связи при его орбитальной 
эксплуатации в существенно нестационарных 
условиях суточного цикла освещенности Солнцем. 
Суммарная мощность тепловыделения приборов 
БА, расположенных на обеих сторонах приборной 
панели - 2 составила 650 Вт. Рассматривалась 
СОТР на основе сети аммиачных нерегулируемых 
НТТ, проложенных в восьми параллельных пло- 
скостях. Эффективные коэффициенты теплопро- 
водности анизотропных элементов нагретой зоны 
приборов БА рассчитывались по тепловым схемам 
элементарной ячейки, представленным в [6] : 
А ѵ =3,3; А.=9,3; Я 1 =9,3 Вт/(м-К). Эффективные ко- 
эффициенты теплопроводности анизотропного со- 
тового заполнителя с ячейкой в виде шестигранни- 
ка рассчитывались по эмпирическим зависимо- 
стям, приведенным в [3]. 

В качестве опорного варианта моделирования 
рассматривались типичные условия орбитальной 
эксплуатации КА в режиме сеанса связи в середине 
срока активного существования (7 лет, Ду=0,26 [3]) 
в точке весеннего равноденствия в суточном цикле 
освещенности Солнцем. Результаты моделирова- 
ния работы ИЛБ в типичных условиях эксплуата- 
ции представлены на рис. 2-4. Характер и диапа- 
зон изменения во времени средних температур об- 
шивок панелей - 1 и 3 позволяет оценить суммар- 
ную эффективность работы СОТР, рис. 2. 

Рост средней температуры обшивок радиатор- 
ных панелей - 3 (0<Э<12, кривые 3, 4) и 1 (12<Т<24, 
кривые 1, 2) более, чем на 10 °С, обусловлен исклю- 
чительно работой НТТ при переносе энергии, вы- 
деляемой приборами БА, и тепла солнечного излу- 
чения на радиаторные панели, находящиеся в ука- 
занные периоды времени в тени. С учетом того, что 
приборная панель - 2 не является радиаторной и 
имеет довольно высокую плотность двухсторонней 
компоновки расположения тепловыделяющих при- 
боров БА, характер изменения средних температур 
обшивок этой панели (кривые 5, 6) также свиде- 
тельствует о том, что в условиях периодической ос- 
вещенности радиаторных панелей Солнцем доми- 
нирующим является механизм переноса тепловой 
энергии с помощью НТТ. Следует отметить, что для 
большинства узлов и блоков радиоэлектронной ап- 
паратуры, используемой в КА, диапазон рабочих 
температур составляет от -50 до +50 °С. 


Т, °С 



Рис. 2. Изменение во времени средней температуры обши- 
вок радиаторных панелей 1 (кривые 1, 2), 3 (кривые 
3, 4) и приборной панели 2 (кривые 5, 6) при орби- 
тальной эксплуатации ИЛБ в суточном цикле осве- 
щенности Солнцем 
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Рис. 3. Изменение температуры приборов БА во времени: 
1) максимальной и 2) минимальной (опорный вари- 
ант); 3) максимальной, режим «Перегрев»; 4) мини- 
мальной, режим «Переохлаждение» 

Известно [8], что увеличение температуры на 
10 °С в области +40.. .+50 °С снижает основные по- 
казатели надежности на 50 %. Поэтому работа РЭА 
в малой окрестности максимальных температур не- 
целесообразна. Анализ динамики изменения су- 
точной максимальной и минимальной температур 
приборов БА в типичных условиях орбитальной 
эксплуатации ИЛБ (рис. 3, кривые 1, 2) показал, 
что могут быть превышены предельные значения 
рабочих температур приборов радиоэлектронной 
БА (+50 °С и более). Кроме того, выявлена доста- 
точно высокая степень неизотермичности (более 
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30 °С) распределения температурных полей по при- 
борам одной панели. Дополнительно об этом мо- 
жет свидетельствовать анализ положения характер- 
ных изотерм наиболее теплонагруженной плоско- 
сти обшивки приборной панели - 2 в момент вре- 
мени достижения максимального значения темпе- 
ратуры приборов, рис. 4. 

С целью оценки теплопереноса в ИЛ Б при эк- 
стремальных условиях его работы, были смодели- 
рованы режимы « Перегрев » и «Переохлаждение». В 
режиме «Перегрев» рассмотрены условия орбиталь- 
ной эксплуатации КА во время сеанса связи в кон- 
це срока активного существования (12-15 лет) в 
точке зимнего солнцестояния в суточном цикле ос- 
вещенности Солнцем. Результаты вычислительно- 
го эксперимента показали, что в подобных усло- 
виях может быть значительно превышен диапазон 
максимально допустимых рабочих температур ра- 
диоэлектронных приборов БА (рис. 3, кривая 3). В 
режиме « Переохлаждение » учтены условия орби- 
тальной эксплуатации КА в начале срока активно- 
го существования в точке летнего солнцестояния в 
суточном цикле освещенности Солнцем. Анализ 
соответствующих результатов моделирования по- 
казал, что в данных условиях приборы БА будут ра- 
ботать без существенного снижения показателей 
надежности (рис. 3, кривая 4). 

Представляет интерес сравнение полученных 
результатов с соответствующими характеристика- 
ми температурных полей, полученными при ис- 
пользовании плоской нестационарной модели [5], 
не учитывающей также совместность протекания 
процессов теплопереноса в различных панелях. Та- 
кое сравнение показывает, что пространственная 


модель приводит к существенно меньшим откло- 
нениям максимальных и минимальных температур 
от средних значений в идентичных условиях тепло- 
вого воздействия как внешней среды, так и вну- 
тренних источников тепловыделения. При этом, 
отклонения составляют для плоской модели до 
25 °С, а для пространственной - не более 12 °С. 
Учитывая вышеупомянутые ограничения, на ди- 
апазон рабочих температур радиоэлектронной ап- 
паратуры КА, можно сделать вывод о существенно 
ограниченных возможностях плоских моделей. 
Следует отметить, что в рассматриваемых условиях 
учет реальных термических сопротивлений пане- 
лей в пространственной модели должен был при- 
водить, на первый взгляд, к обратному эффекту: 
большим по сравнению с плоской моделью пере- 
падам температур. Установленные закономерности 
могут быть объяснены тем, что плоские модели по 
существу не учитывают теплоаккумулирующую 
способность основных элементов конструкции КА 
- панелей. В экстремальных режимах работы те- 
плоемкость конструкции играет значимую роль, 
снижая интенсивность разогрева в режиме « Пере- 
грев » и увеличивая температуру панелей в режиме 
« Переохлаждение » при нестационарном внешнем 
воздействии. Полученные результаты показывают, 
что исследуемые процессы в такой сложной техни- 
ческой системе, как ИЛБ современного КА аде- 
кватно можно описать только при моделировании 
с применением пространственной нестационарной 
математической модели теплопереноса. 

На основании результатов проведенных теоре- 
тических исследований можно сделать вывод о 
том, что характерными для Н-образного ИЛБ со- 






Рис. 4. Положение характерных изотерм наиболее нагретой обшивки приборной панели 2 при достижении максимальной 
температуры приборов 
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временного КА связи являются существенная не- 
стационарность и неоднородность температурных 
полей всех панелей блока. Масштабы отклонений 
предельных значений температур от средних вели- 
чин существенно отличаются от аналогичных дан- 
ных, полученных с применением плоской неста- 
ционарной модели теплопереноса. Следовательно, 
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